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O . Abstract 

o ■ 

O ' We applied QCD Light Cone Sum Rules to estimate power corrections to the 

. helicity-conserving amplitude in the process 7*7 vrvr. We found that above ~ 4 

QhI GeV^ power corrections are numerically small and the twist-2 part dominates. The 

Q^' amplitude can be reliably calculated in this region using models of 2tt distribution 

1^ . amplitudes as an input. We found that the magnitude of the NLO corrections 

'"^ ' depends rather strongly on the normalization of the gluonic distribution amplitude. 
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1 Introduction 

Hadron production in the reaction 7*7 liadron(s) has been a subject of considerable 
interest for a long time, both from the experimental |^ and theoretical p, ^, |^ points 
of view. The key role in QCD description of such processes is played by the QCD fac- 
torization theorem. For example, QCD factorization has been successfully applied to the 
reaction 7*7 ^ 7r° [g g. The F'f'^iQ'^) form-factor data obtained by the CELLO and 
CLEO collaborations are in a good agreement with the available QCD analysis, see for 
example g 0, @. 

Recently it has been proposed 0, 10| to investigate a similar process 7*7 — > tttt when the 



two pion state has a small invariant mass. It has been argued that QCD factorization 



applies to this case as well [|1T| . The resulting amplitude depends on new non-perturbative 
objects, the so-called two-pion distribution amplitudes (27rDA's). They are given by 
matrix elements of twist-2 QCD string operators between vacuum and the two-pion state 



12| . Moreover, 27rDA's can be related by the crossing symmetry to skewed parton 



distributions [0, [Tj] which recently have been subject of considerable interest. 
Furthermore, in the recent paper [|l^ it has been argued that experimental studies of 27r 
production cross-section are possible with existing e~^e~ facilities. 

Formally, dominance of the leading-twist amplitude is guaranteed only at a very large 
For the process 7*7 — > vrvr the bulk of the twist-2 amplitude arises from the handbag 
diagram. In analogy with the 7*7 7r° reaction one expects that the leading-twist 



contribution dominates the amplitude even for moderate values of ~ 4 — 10 GeV^. 
For lower values of , the power-suppressed corrections are certainly important. Note 
that preliminary estimates show that most of the 7*7 mr events which have been seen 
in the CLEO data are in the region ~ 1 — 5 GeV^ 0, |l^. In this region a reliable 
estimate of the amplitude including the power-suppressed contributions is crucial. 
From the theoretical point of view problems encountered in an analysis of power-suppressed 
contributions to two-pion and one-pion production amplitudes are very similar. Terms 
suppressed as 1/Q^ can arise from different space-time configurations. Production of 
states with more than 2 partons by interaction of two electromagnetic currents at small 
transverse distances can be accounted for by the standard Operator Product Expansion 
(OPE) technique. However, as the second photon is real, there is yet another config- 
uration which results in a power-suppressed correction. The real photon can turn into 
hadrons long before the interaction with the virtual one. It occurs at large transverse 
distances between two electromagnetic currents. Such a not-factorizable term is known 
in the literature as the 'end-point' or 'soft' contribution. 

In this paper we use the light-cone sum rules (LCSR) method to evaluate the 7*7 tttt 
amplitude, including the power suppressed contributions. The main advantage of this 
technique is that it allows to take into account both factorizable and non-factorizable 
corrections. Recently, LCSR were successfully applied to describe different pion form- 
factors 0, [H in a Q2 range from 1 to ~ 10 GeV^ 

This paper is organized as follows. In the following section we present the definition of 
the 7*7 — s> TTTT amplitude and set up the notation. In section 2 models of 27r distribution 
amplitudes are introduced. Sections 2 and 3 are devoted to discussion of the LCSR to the 
LO and NLO accuracy, respectively. In section 4 we present numerical analysis. Finally, 
we summarize. The appendix contains definitions of the NLO Wilson coefficient functions. 

2 General definitions 

Kinematics of the reaction 'y*{q)'y{q') ~^ 7T{ki)7T{k2) can conveniently be described in 
terms of a pair of light-like vectors p, z which obey 

p2 = ^2 ^ p . ^ _^ Q 

and define longitudinal directions. Here p ■ z = p^z^. Let P and k denote total and 
relative momenta of the vr meson pair, respectively, 

p2 = + ^2)2 _ p^2^ ^2 _ {ki - k2y = Ami - W^, P ■ k = (2) 
The initial and final states momenta can be decomposed as 



q = p 



2{p-z] 



-z, = -Q^ 



q = — — z, q = 



2{p-z) 



P = q + q' = p + 



2{p-z] 



(3) 



2[p ■ z) 



The longitudinal momentum distribution between pions is described by the variable ^ = 
{k ■ z)/{p ■ z). Alternatively, 

^ = /3 cos dcra , 

where ^cm is the polar angle of the pion momentum in the CM frame with respect to 
the direction of the total momentum P and j3 is the velocity of produced pions in the 
center-of-mass frame 
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The amplitude of hard photo-production of two pions is defined by the following matrix 
element between vacuum and two pions state: 

T^"^ = i j d'^xe-'''-'^{27i{P, k)\TJ''{x/2)r{-x/2)\0) , Q = ^{q - q') (4) 

where J^{x) denotes quark electromagnetic current. Hard photo-production corresponds 
to the limit ^ W"^ > Aq^q where the amplitude @ can be represented as an expan- 
sion in terms of powers of 1/Q. According to the factorization theorem the leading twist 
term in the expansion can be written as a convolution of hard and soft blocks. The coef- 
ficient functions can be calculated from appropriate partonic subprocesses 7* -|- 7 — > g -|- g 

01 Y + 1 9 + g- 

According to the analysis of Ref. [0 to the leading twist accuracy the amplitude T^^ is 
a sum of two terms 

T^\q, q\ P, k) = '-{-gnrTriq. P. k) + T^iq, P. k) (5) 

where {—g^^)T = ( ^^^' ^^p"^^ — gf^'^^ is the metric tensor in the transverse space and k^_l^k'2 
denotes traceless, symmetric tensor product of relative transverse momenta (p7|). 
The leading-order, leading-twist amplitude Tq'^'' describes scattering of two photons with 
equal helicities, related by crossing to the photon helicity-conserving DVCS on a pion. At 
the NLO there is a new contribution T2^'^ from collisions of photons with opposite helic- 
ities, related to photon helicity-fiip contribution to DVCS [^]. In terms of the Operator 
Product Expansion the latter amplitude singles out a twist-2 tensor gluon operator which 
cannot be studied in deep-inelastic scattering (DIS) on a pion (or nucleon) target |T9[] . 
Note that as the amplitude T^'^ is dimensionless, twist-2 amplitude Tq'"'^ depends on 
only logarithmically through the running coupling and QCD evolution effects. To see this 
it is convenient to develop an appropriate power counting in an infinite momentum frame. 
For definiteness we assume that the pion pair is moving in the positive z direction and 
and z~ are the only nonzero component of p and z, respectively. Then the infinite 
momentum frame can be understood as p'^ ^ Q ^ 00 with a fixed {p ■ z) ^ 1. From (^) 
it follows that in this frame {k ■ z) ^ 1 and kj_ ~ Q^. This determines the power counting 
in Q for all twist-2 amplitudes and from we find that Tq'^'" is 0(1) as far as powers of 
Q are concerned. 

In this paper we consider power corrections to the amplitude Tq^'^ only. It is expected to 
be the most important term numerically as it appears already at the Born level. To the 



NLO accuracy one has ||T9[|: 



rp'ynTv rppert 



4tt 



(6) 



Coefficient functions Cg,Cg,Cg can be found in [|T9l. The quark and gluon 27rDA's are 
defined as matrix elements of the hght-cone string operators: 



Nf Y ^0 



(7) 



(7r7r(P, k)\z>'z'' |0) = {p ■ zf / du W^) e'(2«-i)(P-) 



Note that distribution amphtudes VT^) and VT^) depend also on a fac- 

torization scale /X. 

3 Models for 2-pion distribution amplitudes 

In this section we describe briefly the main properties of the distribution amplitudes 
introduced in (|^ and (H) and discuss a model which has been used to obtain estimates 
for the magnitude of power-suppressed corrections. 

Due to the positive C-parity of the pion pair, 2tt distribution amplitudes have the following 
symmetry properties 

^, W^) = -$^(1 - u, ^, W^) = W^), 

(9) 

e, w^) = $^(1 - u, e, W) = -e, w^). 



The factorization scale dependence is governed by the ERBL evolution equations |22 



As it is well known, in the leading logarithmic approximation their solution has a form of 
an expansion in terms of Gegenbauer polynomials: 



(10) 



n=l 
odd 



<I>^K W'\ /.) = 30 u\l - uf Mi. ^) Cl'\2u - 1), 



n=0 
even 



Coefficients Bn and mix under evolution. 

In the next step one expands both distribution amplitudes in the partial waves of the 
final 2-pion system As a result one introduces an expansion of VT^I /i) and 

y4„(^, W'^l /i) in terms of Legendre polynomials Pi{C) [[ 
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n+l 



B^{i,W^\^^)= Y.B^l{W^\^^)m) 



1=0 
even 
n+2 



(12) 



A4^,w'\fi) = J2A^iiw'\fi)m) 



1=0 
even 



Additional constraints on 27r distribution amplitudes are provided by soft pion theorems, 
see 



12 



^ = 1,W^ = 0) = ^^{u, e = -1, W^ = 0) = 0. 
^ = 1,W^ = 0) = ^^{u, ^ = -1, = 0)= 0. 



(13) 



Finally, crossing symmetry allows to relate moments of distribution amplitudes to forward 
matrix elements of twist-2 operators which determine moments of pion quark and gluon 
structure functions, see [jl2| for details. One finds, e.g. 

du{2u - 1) ^ = 1, iy2 _ 0) = --^M«(i - e^) 

Jo Nf 

du ^ = 1,W'' = 0) = --M^{1 - 
Jo 2 

where M*^ and M*^ are momentum fractions carried by quarks and gluons in a pion: 

M^{^)=[ duuy2iQn{u, n) + q^riu, fx)) 
Jo „ 



(14) 



M'^{fi)= I duugT,{u, fi) . 



(15) 



At asymptotically large — > oo only the lowest terms in (|l3), (0) contribute. Combining 
constrains ([T^ with ( [T^ ) one easily finds: 



$^,(n, C, W' = 0) = -Ibu'il - «)^M,^,(1 - e) , 

(w, C, W' = 0) = -30«(1 - «) (2« - 1)^M«(1 - a 



(16) 



where 



AT, 



Ar/ + 4CF' 



m: 



G 



4Cf 



AT/ + 4Cf 



(17) 



For W"^ 7^ one writes 



2 ;\/rG 



$^,(«, e, W) = -301.(1 - tz) {2u - 1)— M«i?(e, W^') 

7V/ 



;i8) 



Function B{^, W ) is related to coefficients Bi and Aq from ([121) 



3 A^f 



5i(^,PF'|/i' = oo 



(19) 



In the limit ^ one finds from (|TD that B{^, = 0) = {1 - ^ 



In numerical calculations presented in this paper we have used a model which retains sim- 
ple analytical form of asymptotic distributions amplitudes ([l^), but incorporates nontriv- 
ial information about pion structure at a scale of the order of a few GeV^ |15]. Assuming 
dominance of the lowest conformal wave one finds 

1 (20) 

<l>^{u,^,W^) = -30m(1-m) (2m-1)— M^(/i2)5(e,Vr). 

Nf 

As at present a little is known about 2tt DA's, the main motivation beyond the model 
(pop is that it incorporates all constraints arising from crossing symmetry and soft pion 
theorems and has a simple form which makes its treatment in numerical calculations 
easy. Note that dominance of the lowest conformal wave seems to be phenomenologically 
justified in the case of the single pion DA. 

This model differs from the asymptotic DA (|1^) only in values of momentum fractions 
M*^(/i^) and M'^(/i^). Their scale-dependence is given by 

MW) = M« (l + L(Q2)i?(/.2)) , R{t,^) = ^''^^2Q^^' ^ (21) 

IVla' 



where L is the usual evolution factor: 



Obviously, that M^{Q^) = 1 - M^iQ^). 

Explicit expression for -B(^, W) can be obtained using the Watson theorem |]12|. In cal- 
culations considered here B{^,W) enters as a Q^-independent factor and therefore its 
explicit functional form is not important for considerations of power-suppressed correc- 
tions to the amplitude. To remove this factor from numerical calculations we will consider 
the Q2_dependence of the ratio T(?'^^(Q^^, W^)/T^', where 

= TriQ' = oo, ^, =Y.<t = E W') . (23) 

is asymptotic value of the amplitude. 

In Fig. in we compare this ratio for the NLO amplitude T^^''\Q'^,^,W^) for different 
models of 27i distribution amplitudes: asymptotic (|r^) and its minimal extension (PP]). 
Note that there is an ambiguity due to the scale-dependence of the NLO amplitude. In 
order to estimate this uncertainty we make the following choice for the scale /x: 

^'^ = kQ"^ + , (24) 

with = 1 GeV^ and a parameter k which will be varied in the interval 1/5 < k < 1. 
Such a choice is motivated by an observation [|l8l that a typical virtuality of a propagator 



in a perturbative, exclusive amplitude is given by a weighted sum of the hard scale 
and the infrared cut-off. In the LCSR approach the latter role is played by the Borel 
mass M^, see the next section. Evaluating the NLO amplitude (0) we have neglected 



numerically small NLO corrections |23, 24, 25 1 to evolution of distribution amplitudes in 



the MS scheme. 



One finds that in the model based on asymptotic distribution amphtudes the NLO cor- 
rection is numerically much larger than in the second model, where amplitudes have the 
asymptotic form but different, scale-dependent normalization. It can be understood by 
observing that the gluon distribution amplitude gives the main contribution to the NLO 
correction. At of order of few GeV^ its normalization in (|20|) , based on the momentum 



fraction carried by gluons in a pion according to the GRV parametrization [26] is much 
smaller than the asymptotic one. Strong sensitivity of the NLO corrections to gluon dis- 
tribution amplitude is an interesting feature of the process considered here. As it follows. 



we have to revise one of our conclusions from |]19[, about the size of the NLO correction. 
Contrary to our previous claim, its magnitude turns out to be rather model dependent, 
and therefore it is difficult to make a trustworthy prediction unless distribution amplitudes 
which enter the NLO correction, in particular the gluon one, are sufficiently constrained. 

4 Light-cone sum rules method: the LO approximation 

In this section we discuss the Light Cone Sum Rule (LCSR) for the amplitude Tq'^'^ . Our 
procedure closely follows investigation of the photon-pion transition form-factor F^t^{Q^) 
inRef. 0. 

The first step to obtain a LCSR for the amplitude with one real photon T^^'^ is to consider 
an amplitude Tq ^'^ where both photons are off-shell and have large virtualities: 

7*{q) + i\q) - MP. k), Q\ Q'^ » a^^d (25) 

and find its dispersion representation in the variable q'^. 

In the general kinematics momenta q and q' can be represented in terms of vectors p and 
z as 

g = -TTT^^P + <rz,q' = -Q' q' = "^^^^ P + az,q" = -Q'\ (26) 
2a[p ■ z) 2a[p ■ z) 

Coefficients a and a are functions of kinematical invariants. Using momentum conserva- 
tion one finds 

— w/^ _ 

^ ~ 2(p-z) 



a 



1 (27) 

-{q . {q - q') -Vx), X = {q ■ q'f - q\" . 

2{p ■ z) 



With the help of the factorization theorem one can write '^'^ as a convolution of hard 
and soft blocks. Virtuality g'^ enters the hard part only, where we neglect as compared 
with Q2 and Q'^. 

With g'^ 7^ the two-photon amplitude admits naturally a reacher Lorentz structure 
than the original amplitude with one photon on-shell. In addition, by splitting T^'^ into 
Lorentz tensors and invariant coefficient functions it is advantageous to avoid kinematical 
constraints for the latter. Constraints imposed on coefficient functions result in constraints 
on the form of their dispersion representation. 



To this end we rewrite the transverse metric tensor in terms of momenta q and q' [Q: 

R'^iq. q) = i-gnr = -g'" + ^{q- q'{qW' + q'^q") - q'q'^q"' - q'^'q"). (28) 

with X defined in (0), and introduce a variable uj 

2P-q _W^ g2 _ ^ g2 _ 
g2 Q2 q2 q 



With these definitions it is convenient to introduce Tq as 

T'-" = ^Y^'^iq, q'W^^iQ. q\ i.W) + . . . (30) 

where elhpsis denotes other possible Lorentz structures. Now, in the limit —q'^ 
[uj 1) Tq defined here, goes smoothly into Tq'"'^ given by (^. On the other hand, 
factor oj"^ cancels singularity present in the tensor R^^" as — g'^, W"^ — » (cu — 0) due 
to the 1/X term, hence no constraints have to be imposed on '^'^ . 
To the LO accuracy and keeping W"^ = in the hard block the amplitude Tg '"'^{Q, q' , C,, W) 
is the only one which appears in T^'^ and one obtains: 



q Jo 1- 

q Jo (1 — x)Q^ — xq'^ 



XUJ 



Note that the above expression depends on W through the 2-pion distribution amplitude 
(0)- 

The LO result can be easily converted into a dispersion integral over g'^ with s = xQ'^/x 
being the mass of the intermediate state. Using symmetry properties of the quark DA 
one finds 

To'*-'^(Q, g', W) = r (32) 
JO s — q'^ 

where 

PoiQ,s,^,W) = 2iJ2el) f dx5{x--^)x''^'^{x,i,W^) (33) 

The next step is to rewrite the dispersion relation in the g'^ channel assuming that the 
spectral density can be approximated by contributions of low-lying hadron states p, uj and 
a continuum of higher-mass states with an effective threshold Sq: 

(g,g,e,W^) = V2/, +/ ds -— ^ . (34) 

171 p — q J so s — q 

Here we have introduced the following notation for the matrix elements: 



(35) 



m„ 2 



e*, ey are the polarization vectors of the p meson. To include tu-meson we adopt the 
following approximate relations 



rriu, ~ m 



pi 



3/. ^ /p, Tr{Qi W) ^ 3Tr{Q, W) (36) 



which follow from the quark content of p and uj and from the isospin symmetry. 
The right-hand-side (RHS) of equation (|3^ ) involves two unknown functions: the form- 
factor Tq^^ and the spectral density Assuming quark-hadron duality one can esti- 
mate the continuum spectral density as 

p^°"*(Q, s, W) = 9{s > so)po{Q, s, W) (37) 



where pq{Q, s,^,W) is the spectral density given in (|33|). By keeping —q'^ large and 
combining (^Tf), ( ^4]) and (|3lD one obtains a LCSR for the form-factor Tq^^: 

mj^ — q'^ Jo s — g'^ 

After perfoming Borel transformation in —q'^ one finds (M is the Borel mass): 

V2f,TriQ, W) = r dspoiQ, s, e, l^)eK-)/A^^ (39) 

Jo 

Finally, substituting (^7[) and ( p9| ) into (0) results in 

To'^*^^(Q, W^) = r ci.po(Q, e, Ty)eK-^)/-^ + r ^40) 

ml — 0^ JO Jsn s — a 



The above representation allows to perform an analitical continuation to the point g'^ = 0. 
In this way one arrives at a LCSR for the amplitude with one real photon: 

TriQ, ^,W) = ^ r dspoiQ, s, iy)eK-)/M^ + r ^po(Q, s, ^, W) (41) 
Jo J so s 

It is convenient to rewrite this formula going back to an integral over fraction x = Q^/ {s + 
Q"^). Introducing xq = Q^/(so + Q^) and x = 1 — x one finds 

TriQ,^,W) = g(E 4) £ dx^^ix, W^) exp I I + (42) 



"p q ^0 



„ ^ JO 1 — X 



Consider now the limit oo. In this limit xo = 1 — So/Q^ + 0{1/Q'^). The integration 

region in the first term in the RHS (|42|) shrinks to the point x = 1 and one obtains 



dx^^ix, IV^) exp I ^^f I = #^$g(l, IV^) tdxxe'^+Oiimm 



rn?p Jxo ' ' [ xM^ J Q^m^ ^ ' ' 7o 

where $^(l,e,W^') = ;^$«(x, l^)U=i. 



As it has been discussed at length in the hterature, equation P3| ) can be interpreted 
as the so-called "end-point" contribution which arises from large transverse dinstances 
between two photons in the hard block [18, 27|. In general, the 27r distribution amplitiude 
ly^) depends also on the factorization scale fi which separates large and short 
distances. As x ~ 1 the quark virtuality (the magnitude of the quark denominator in 
equation (pi])) is of order of the Borel mass M. As it follows, the two pion DA in (|43|) 
should be evaluated at a low normalization point, of order of M. 

In Fig.^ we show an average value of the momentum fraction x in the integral (|l2|) 
calculated as a function of Q^. One observes that the mean value of x in (^) is indeed 
close to 1. 

As — > c>o the second term in the RHS of (^21) gives 



dx^:!—^'^{x,i,W^) 

1 — X 







1 

dx-^^'^{x,i,W^) + 0{l/Q'') 

\ — X 



(44) 



It reproduces the leading order, leading twist factorization formula when g'^ = in (|31|) 
and provides correct asymptotics for very large Q^. The power correction is suppressed 
as 1/Q2. 

The LO sum rule (^21) results in an expression for the amplitude T^'^'^ which includes 
contributions from both the "end-point" region a; ~ 1, associated with large transverse 
distances, and from small transverse distances where the qq pair is created by two photons 
in a compact configuration. From the sum rule it follows that the "end-point" contribution 
is suppressed by as compared with the LO factorization result, i.e. it has the same 
order as factorizable higher twist corrections. Despite formal power suppression, the "end- 
point" contribution can be numerically important for realistic values of Q^. Our numerical 



analysis suggests that the sum rule ( ^2|) can be applied for description of the amplitude 
T'^'^'" starting from moderate momentum transfers > IGeV^. 



Radiative corrections 



In principle, the sum rule ( ^2|) can be improved in a twofold way, by including the NLO 
tts and higher-twist corrections to the spectral density (|33D. In this section we consider 
the NLO contribution. Taking into account higher-twists requires knowledge of the cor- 
responding 2- pion distribution amplitudes which are not known yet. 



One-loop corrections for the real photon case have been considered in In current 

situation one should calculate the coefficient functions in kinematics when both photons 
are vitual. One should also keep in mind that at the NLO diagrams with gluons enter the 



game. As in [TU| one can use crossing symmetry to derive coefficient functions from the 
corresponding coefficient functions in the DVCS kinematics, as computed in |2^, |2^, |3^ . 



The calculation is straightforward. The NLO amplitude can be written in the standard 
form: 



(45) 



T^*''^iQ,q',^,W) = Y.el [I dx^Qix,^,W')[c'^ix,uj) + ^Clix^u^f^) 

dx^''ix,^,W') \^Cliz,co,fx) 



Coefficient functions C° Cl and are collected in the Appendix. 



As in the LO case, virtuality q'^ enters only through the variable uj. It is convenient to 
rewrite equation (^51) in the form which resembles the structure of the LO answer (|3TD : 







1 — XUJ 



with 



p{x,i,W') = p,{x,i,W') + 



2^ , 



An 



(46) 



(47) 



As in equation (|3T| ) 

poix, e, W) = 2(^ ej)<l>«(x, W) . (48) 
The NLO correction pi{x, ^, PF^) can be obtained from the corresponding coefficient func- 



tions ([451) : 

^ el /' rfx [<l>«(x, e, W^')-Im,,2C^i(x, a;) - $^(x, ^, W')-lm^,2Cl{x, u) 

— JO L TT TT , 

= mV(m,^,W^^)I«=i/.. 



(49) 



where Imq/2 denotes the imaginary part with respect to the variable q , considered here 
as a complex variable with positive both real and (infinitesimal) imaginary parts. 
With such a definition the structure of the NLO LCSR for the amplitude Tq'^^ remains 
the same as the LO one: 



rriQ, ^,W) = % f dxp{x, e, W) exp 



xm'p — xQ^ 



xM^ 



+ 



dx- 



X 



1 — X 



p(x,e,i^').(5o) 



As in the case of the LO sum rule (^21), as — *• oo equation (^) reproduces the pertur- 
bative expansion of the Tq^'^ amplitude, including the NLO corrections The power 
correction is suppressed as 

As an illustration, we quote now the form of pi(x,^, W^^) in the case of the 'minimal 
model' (BOl) 



pi(x,e, W^') = Y.4B{i,W^) {ln(/iVg')pio(x,/i') +Pii(x,/i2)} ^ 
pio(x, p2) = (-40)x(l - x){2x - l)i?(p^), 

Pll(x,p^) = Pas{x) + R{p^)[pas{x) + D{x)] . 

Where -R(p^) is defined in (|21|). Function D{x) is a shorthand notation for: 

D{x) = - — [2x(l - 6x) + xx{2x - 1)(31 + 121n(x/x))] 
3 

Taking for the asymptotic value (|1^) results in _R(p^ = oo) = and 

p,{x,^,W')=Bi^,W')pasix) 



(51) 



(52) 



(53) 



where 



Pasix) = ^^^-^^(1 - ^)(2a: - - 17 - 3 ln^(x/x)) (54) 
In this case pi has no dependence on /i. 

Before evaluating the sum rule one has to provide an estimate of the factorisation scale /i. 
This is a standard problem in a calculation based on the fixed-order perturbation theory. 
Note that the scale of perturbative expansion /i^ ~ is different from the characteristic 
scale of soft, 'end-point' contributions /i^ ~ M^. As a consequence, the correct treatment 
of the sum rule requires applying various normalization scales to various terms in (|50D. 
To avoid confusion, we propose the following procedure. One separates the perturbative 
contribution and writes the final formula in the form: 

TriQ, w) = Tr\Q, e, w) + rr-'^'^XQ, w) . (55) 

Tq"^^* is the perturbative amplitude (|^), which can be represented as: 

Tr\Q.^,W) = f'dx-^p{x,C,W'). (56) 

JO 1 — X 

The scale in this term is set by the hard photon virtuality Q^. The second term should, 
as a matter of fact, be considered as the proper light-cone sum rule result: 

JXQ 

Here the integration region is restricted to xq < x < 1. As discussed in the previous 
section, this term should be evaluated at a low normalization point ~ M^. 
Of course, when evaluated at the same normalization scale, sum of (|56|) and (|57| ) repro- 
duces the original sum rule (|50D. 




xm. 



xQ' 



X 



a;M2 



1 — X 



(57) 
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Numerical results 



In the subsequent numerical analysis the following input has been used: the threshold 
1.5 GeV^ has been taken from the two-point sum rule 

0.5 



parameter sq 



This sum 
0.8 GeV^. In the 



rule is reliable for the corresponding Borel parameter M|pj 
light-cone sum rules it should be larger to compensate for the fact that the effective 
expansion parameter is given by the inverse power of xM"^. In this case a typical choice 
is ~ M^pJ{x) |l3, HI- We assume 0.6 < < 1.2 GeV^ as a reasonable interval. In 
addition we have checked that changing sq by ±0.2 GeV^ does not produce any sizeable 
effect. 

We remind the reader that in the present investigation we have neglected higher twist 
contributions to the sum rule for Tq '^^ . They usual, suppressed by additional 

powers of the Borel parameter. To obtain their contribution one should know two-pion 
distribution amplitudes of higher twists which have not yet been studied. In the case 
of photon-pion transition form-factor [|17|, the contribution of such terms to the non- 
perturbative (power suppressed) part is about 30 — 35% of the leading twist contribution 



for the low ~ 1 — 2 GeV^. Qualitatively, the same picture should hold also in the 
present case. 

All numerical results have been obtained with the model ( 120|) for 2tt distribution ampli- 
tude. We have used the quark momentum fraction M^{0.6 GeV^) = 0.63 taken from the 
GRV parametrization [Q. We have kept Nj = 4 which results in = 0.43. 
In all calculations we have used Aqcd = 0.204 GeV and one-loop running coupling 
asilGeV^) = 0.47. 

In the numerical analysis one has to specify the factorization scale. Since after subtraction 
of the perturbative amplitude Tq°^~'^^^^{Q,^,W) is dominated by the "end-point", soft 
contributions, we apply here a fixed, low scale /i^ ~ M^. 

To estimate the ambiguity due to the scale dependence in Tq^'^^{Q,^, W) we parametrize 
/i^ according to: 

= kQ^ + m2 , (58) 
with a parameter k which varies in the interval 1/5 < k < 1. 

We consider values of in the physically interesting interval 1 GeV^ < < 10 GeV^. 
To cancel the influence of the overall factor W^), in the following we always plot the 
ratio of Tq''"^(Q^,^, W'^) to its asymptotic value at = oo (p3|). 

The Borel parameter dependence of the ratio j'°°'^~P*'''*/j'as jg gj^Q-^^ pjg |^ find 
that the dependence is sufficiently fiat to justify the use of LCSR, although for lower 
values of Q^, where power corrections are increasingly important, dependence on the 
Borel parameter becomes stronger. Assuming that the uncertainty arising from the Borel 
parameter dependence should not exceed 30 % for the j'°°'^~p'=''*^ estimate that the sum 
rule ( ^21) provides a reasonable description of Tq^^ starting from > 1 GeV^. 
In Fig. § we show the ratio Tq"^'"^"^^^ (Q'^) /T^ obtained from the light-cone sum rule (|57|) 



for different values of the Borel parameter M^. One observes that in the whole interval of 
considered here the sum rule calculation is rather stable with respect to variation of the 
Borel parameter within a reasonable interval. In accordance with expectations, the non- 
perturbative correction becomes smaller with increasing Q^. By fitting a simple formula 
we found that for the value of the Borel parameter in the middle of the interval, = 0.9 
GeV^, and in the region 1 < < lOGeV^, jj^°°~P''''*y'j'as ^^^^ parameterized as 

™-pert.^2w™s -l-5 + 0.05Q^ + 0.015g^ 
with Q"^ in units of GeV^, see the solid line in Fig. ^. 

In Fig. 1^ we present perturbative and non-perturbative contributions to the ratio Tq'^'^ /T^ 
as a function of Q^. j'^°'^~v<^^'^ jg calculated with = 0.9 GeV^. The non-perturbative 
corrections is numerically significant only in the region < 4 GeV^. For higher values 
of the amplitude is dominated by the NLO leading-twist contribution. 



7 Summary and conclusions 



The main result of this paper is the numerical estimate of power-suppressed correction 
to the leading-twist helicity-conserving amplitude of the process 7*7 tttt. Light-cone 



sum rules technique allows to circumvent difficulties due to non-factorizability of the 
power-suppressed terms. Although formally our analysis is not complete, as we have 
neglected contribution of higher twist operators to the amplitude of two-pion production 
in a collision of two virtual photons, we believe that the general picture is reliable, at 
least qualitatively. Power corrections are increasingly important with decreasing for 

< 4 GeV^, and become about 50% of the leading- twist amplitude at = 1 GeV^. 
Our final result for the helicity-conserving 7*7 tits: amplitude is shown in Fig. The 
grey bound indicates uncertainty due to the linearly combined Borel parameter and the 
factorization scale variations. One finds that starting from around 4 GeV^ the twist-2 
contribution approximately saturates the amplitude. This observation suggests that the 
cross-section of the process 7*7 tttt can be accurately predicted in QCD for given 
models of 2ti distribution amplitudes. 

Assuming dominance of the lowest conformal wave, we have found that helicity conserving 
amplitude is very sensitive to the normalization of gluonic 27r distribution amplitude. This 
observation, combined with crossing, makes it plausible to use 7*7 — >• tttt to constrain 
momentum fraction carried by gluons in a pion. 
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8 APPENDIX: Coefficient functions 

To the NLO accuracy the amplitude can be represented as a standard convolution 

of 27r distribution amplitudes with the hard scattering coefficients: 



(A.l) 



where u has been defined in (^). The coefficient functions in the MS scheme read 

2x 



C'(x,u) 



1 — xu' 



(A.2) 



Cl{x,uj,fi'/Q') 



Cl'{x,u;) 



(1 — xuj)uj 



ln(l -cj; 



2(1 -cj) 
xu'^ 



1 1 



1 — xu u 



+ 



1 — XUJ 



XUJ^ 



XXU)^ 



XUl 



X 



UJ 



1 — XLJ 



(A.3) 



[1 — xuj)uj 



2 3 3 
-9 H ln^(l -xu) ln(l - xcu) H ln(l - lu) 

XUl XUl XUl 



-— In'ri-cu) 

XUJ 



3 

+ 3-^ln(l - XLv) 



XXUJ^ 



■ln(l - a;) + 



+ (x - 1 - l/u) — ^ ^ + (1 - X + x/o;) — _ ^ ^ 



(A.4) 



(A.5) 



+ (1 - a; + [1 - xujf) ln(l - xu;) - (2 - ujx^ - ujx'^){l - uj) ln(l - , 



C]^(a;,u;) = ,/ , |ai(x, a;) ln(l - a;) + a2(x, a;) ln^(l - (A.6) 



ci(x, a;) ln(l — xuj) + ci(x, a;) ln(l — xuj) + 



€2(2;, a;) ln^(l — xui) + €2(2;, a;) ln^(l — a;a;)| , 



ai{x,uj) = S + A(jj{x-3-x'^)+Auj'^{l-x + x'^), 



a2{x,uj) = -2 + uj{3-2x + 2x'^) -u^{l-2x + 2x'^), 



ci{x,uj) = -8 + 4uj{l + 2x) -2uj'^x{l + x), 



C2(x,oj) = 2-a;(l + 2x) +a;V. 

Here we used the shorthand notation x = 1 — x. Note that physical amphtude does not 
have a singularity (pole) when a; — > and therefore all coefficient functions must be well 
defined in this limit: 



C}{x,u;,i,yQ')^CF 



\n{pLyQ')-{2x-l) + {l-2x) 



+ 0{u;), 



(A.7) 
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Figure 1: NLO results for the Tq^'^^{Q'^)/Tq^ as a function of Q^. The upper plot corresponds 
to the mimimal model with M*5(lGeV^) = 0.60 and the lower one to the asymptotic choice 
M*5 = 0.43. The scale is taken as /x^ = kQ"^ + M^, = 1 GeV^. Gray bands show variation 
of the NLO prediction when k is varied between 1/5 and 1. Solid lines correspond to k = 2/5. 
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Figure 2: Average momentum fraction function Q^. We take 

So = 1.5 GeV^, see explanation in the main text. 
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Figure 3: Ratio T^^'^'-^'^'yT^' as a function of for = i GeV^ (short dashed), = 3 QeV^ 
(long dashed) and Q"^ = lOGeV^ (solid line). 
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Figure 4: Ratio TjJ^""^ pert^jTas ^ function of Q^. Short- and long-dashed lines correspond to 
= 0.6 and 1.2 GeV^, with /x^ = M^, respectively. Solid line (M^ = 0.9 GeV^) represents 
parametrization (^). 
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Figure 5: The ratio Tq^'^ /Tq^ as a function of (solid line). Short-dashed line represents 
T^^^/T^ with fi^ = 2/5 Q2 _^ iGeV^. Long-dashed hue represents t^^'^^'^^''^ /T^ with /i^ = 
= 0.9 GeV^. Solid line is the sum of both. 



1 . 1 




. 5 



2 4 6 8 10 GeV 2 

Figure 6: The LCSR prediction for the ratio Tq'"'^ /T^^ as a function of Q^. The grey band 
shows the sensitivity of our result to variation of the Borel parameter within 0.6 < < 1.2 
GeV^ and factorization scale according to formula (|5^). 
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Figure 1: The LCSR prediction for the ratio T^'"'^ /Tq^ as a function of Q"^. The grey band shows the 
sensitivity of our result to variation of the Borel parameter within 0.6 < < 1.2 GeV^ and factorization 
scale according to formula (??). 
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